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Abstract 

We investigate the effect of a nondegenerate quadratic nonlinear dimeric 

impurity on the formation of stationary localized states in one dimensional 

systems. We also consider the formation of stationary localized states in a 

fully nonlinear system where alternative sites have same nonlinear strengths. 

Appropriate ansatzs have been chosen for all the cases which reproduce known 

results for special cases. Connection of the stability of a state with its energy 

is presented graphically. 
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I. INTRODUCTION 



The discrete nonlinear Schrodinger equation which describes a number of phenomena in 
condensed matter physics, non-linear optics and other fields of physics [1-6] may be generally 
written as 

dC 

^~dt~ = ~^ m ^ m ^ ^ m \)Cm + Kn,m+lCm+l + Kn,m-lCm-l 

where V m , m+l = V^ +l m ; and m = 1, 2, 3, ...n. (1) 

In eq.(l) Xm is the nonlinearity parameter associated with m-th grid point and Kn,m+i is 
the nearest neighbor hopping matrix element . Since, J2 m \C m \ 2 is made unity by choosing 
appropriate initial conditions, \C m \ 2 can be interpreted as the probability of finding a particle 
at the m-th grid point. The nonlinearity term f m (\ C m |) arises due to the coupling of the 
vibration of masses at the lattice points to the motion of a quasi-particle in the high frequency 
lattice vibration limit [3]. The DNLSE in other physical context is known as the discrete 
self trapping equation (DST). A number of works on the DNLSE and the DST have been 
reported [7-10]. 

As for the application of the DNLSE particularly in condensed matter physics, we cite 
among others the exciton propagation in Holstein molecular crystal chain [2]. In general, 
the exciton propagation in quasi one dimensional systems [11] having short range electron 
phonon interaction can be adequately modeled by the DNLSE. Other examples include non- 
linear optical responses in superlattices formed by dielectric or magnetic slabs [12] and the 
mean field theory of a periodic array of twinning planes in the high T c superconductors [13]. 
We also note that the vibration in the nonlinear Klein-Gordon chain in certain restrictive 
conditions can be described by the DNLSE [14]. 

One important feature of the DNLSE is that this can yield stationary localized ( SL ) 
states. These SL states might play a relevant role in the nonlinear DNA dynamics [15] and 
also in the energy localization in nonlinear lattices [16]. It has been shown that the presence 
of a nonlinear impurity can produce SL states in one, two and three dimension [17-24]. The 



2 



formation of SL states due to the presence of a single and a degenerate dimeric nonlinear 
impurity in a few linear hosts has been studied in details [20]. The same problem is also 
studied starting from an appropriate Hamiltonian [23]. The fixed point of the Hamiltonian 
[21-24] which generates the appropriate DNLSE can also produce the correct equations 
governing the formation of SL states. We further note that the appropriate ansatz for the 
dimer problem has been obtained in our earlier analysis [20]. Furthermore, the formation of 
intersite peaked and dipped stationary localized states has been studied using the dimeric 
ansatz [24] . The formation of SL states in a perfect nonlinear chain containing one nonlinear 
impurity as well as a degenerate nonlinear dimeric impurity has been studied [24] . In another 
study the formation of SL states in a perfect nonlinear Cayley tree and in a linear Cayley 
tree with dimeric impurity has been considered. 

So far degenerate nonlinear dimer impurity has been considered . Naturally question 
arises what happens if the impurities in the dimer are nondegenerate . We therefore plan to 
study the effect of nondegeneracy in impurities on the formation of SL states. Three linear 
hosts are studied, namely, a perfect one dimensional chain , a Cayley tree and a linear chain 
with a bond defect . We also consider the formation of SL states in a fully nonlinear chain 
having a regular binary alloy composition. Appropriate ansatzs have been fitted in all the 
cases to obtain the results. 

The organization of the paper is as follows. In sec. II we consider the effect of nonlinear 
nondegenerate dimer on the formation of SL states in linear systems. Sec. Ill deals with 
the formation of SL states in a fully nonlinear system and the stability analysis of the SL 
states. Finally we summarize our findings in Sec. IV. 

II. NONDEGENERATE NONLINEAR DIMER IMPURITY 

We consider a system of one dimensional chain consisting of a nondegenerate nonlinear 
dimer impurity of the kind x|C|°" (power law impurity) where a is arbitrary. The dimeric 
impurity is placed at sites and 1 of the chain. The dimer is nondegenerate in the sense 
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that nonlinear strengths at sites and 1 are different and denoted by xo and xi- Therefore, 
the relevant Hamiltonian for this system is 

1 CO y 1 

h = - 9 p m c + i + ^.] + ^[c'oC'r+^ + ^|c'or +2 + ^-ic' 1 r +2 . ( 2 ) 

2 m =-oo 2 a + 2 a + 2 

Hopping matrix element between sites and 1 is V and others are taken to be unity. The 
Hamiltonian with V = 1 describes a perfect chain, V < 1 describes a system of perfect 
Cayley tree while V > 1 describes 1- dimensional chain with a bond defect between sites 
and 1. Furthermore, all the systems contain a nondegenerate nonlinear dimeric impurity 
spanning the zeroth and first sites. In case of the Cayley tree V takes the form of where 
K stands for the connectivity of the Cayley tree. The relevant transformation of the Cayley 
tree into an effective 1-dimensional system has been already presented [24]. Since we are 
interested in the possible solutions for SL states due to the presence of dimeric impurity we 
assume 

C m = cf> m e- iEt (3) 
The presence of a dimeric impurity in the system suggests us to consider the following form 

of (j) m 

<P m = [SgniE^r- 1 ^ m > 1 

Hm , = [^n(£)??]H0 o m < (4) 

We further note that the above form of m can be derived from Green's function analysis 
[20]. Here r] e [0, 1] is given by rj = \[\E\-^E 2 -A}. We further define (5 = o /0i G [-1, 1]. 
Now normalization condition, J2 m=-oo I Cm 1 2 — 1 together with eqs.(3) and (4) give 

1 - n 2 

Using eqs.(3), (4) and (5) in the Hamiltonian H, we obtain the effective Hamiltonian of 
the reduced dynamical system given by 

CT+2 

H eff = Sgn(E) V + Vp\=^ + _L_ (^^j ' (x(l + \W +2 ) + 5(1 - |/T +2 )) (6) 
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Here we have defined x an d 5 as x = (xo + Xi)/2 and S = (xo — Xi)/2. The Hamiltonian 
consists of two variables namely (3 and r\. SL states can be obtained from the fixed point 
solutions [21-23] of the Hamiltonian, H e ff, of the reduced dynamical system. The fixed 
point solutions satisfy equations dH e ff/d(3 = and dH e ff/dr/ = 0. 
Equation dH e ff/d(3 = gives 

' " /»[(* + 6) - W(x -S)\ 
Equation dH e ff/drj = gives 

„ = Sgn(E)/3[( X + S)-\l3nx-S)] (R) 
V V[(x + 5)-\(3\^( X -S)] ' { ' 

For fixed values of x an d ^ number of solutions satisfying eqs.(7) and (8) simultaneously 

will give the number of possible SL states. This can be obtained for arbitrary a but we will 

consider a = and a = 2. The reason behind taking a = case is that the result is already 

known and hence the appropriateness of the ansatz can be verified, a = 2 is considered 

because it is physically more relevant. 

1. (7 = 

For <T = 0, the system reduces to a perfectly linear 1-dimensional chain with two static 
impurities xo an d Xi & t sites and 1 respectively. In this case eq.(7) directly gives 
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In passing we note that , when 5 = i.e. the dimer impurity becomes linear and degenerate, 
(3 — ±1 are only permissible solutions. This is consistent with our earlier work [20]. On 
the other hand 5^0 necessitates f3 ^ ±1 for permissible solutions. Substituting eq.(9) in 
eq.(8) we obtain 

(10) 



X Sgn{E)TvVS 2 + V 2 ' 
Eq.(10) along with signature of E gives four expressions for 1/x- From each expression we 
will get a critical line in the (%, S) plane. The critical lines can be obtained by putting r] — 1 
in eq.(10). The equations describing the critical lines are 



(1) 



1 + V5 2 + V 2 




i-Vs 2 + v 2 




i - Vs 2 + v 2 



-i + Vs 2 + v 2 . 



(11) 



These critical lines are shown in fig.(l) for V = 1, in fig. (2) for V = 0.5 (which corresponds 
to a Cayley tree with K = 4) and in fig. (3) for V = y/2. There are several regions in the 
(x, 5) plane bounded by the critical lines in all the figures. There are regions containing no 
SL state, one SL state appearing above the host band, one SL state appearing below the 
band, two SL states with one below and one above the band and two SL states appearing 
above as well as below the band. Since our results for V — 1 and V = 0.5 agree with known 
results [25], the validity of our starting ansatz is established. 



We now investigate the case where nonlinearity power a = 2 . Physically a — 2 can 
appear when the masses at the lattice points are treated as Einstein oscillators. From eq.(7) 
we again note that for 5 — 0, (5 — ±1 are permissible solutions. On the other hand, for 
5 7^ 0, f3 = ±1 are not permissible solutions. From eqs.(7) and (8) we obtain the relevant 
equation for SL states. The equation is given as 



where a = 5/x- Using eq.(12) we have obtained the full phase diagram of SL states in 
(x, 8) plane for V — 1 and this is shown in fig. (4). There are various regions in the phase 
diagram numbered according to number of possible SL states in the region. Along the 5 = 
line, we note that for x > 0, there are three critical values of x, namely, 1, 8/3 and 8. These 
critical values are in agreement with the values obtained in our earlier work [23]. In this 
case, maximum number of SL states appearing is five. With the use of same eq.(12), we 
have also obtained the phase diagram for SL states in (x, 5) plane for a Cayley tree with 
connectivity K — 4 (V — 0.5). This is shown in fig. (5). Maximum number of possible SL 



2. a = 2 



1 _ g[(l + a) - (1 - a)(3 2 } [V 2 [(l + a) - (1 - a)^] 2 - (3 2 j(l + a) - (1 - a)(3 2 } 2 } 
X ^[l-/3 4 ][(l + a0-(l-a0/3 4 ] 2 



(12) 
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states in this case is found to be seven. The critical values of x along 6 = line are again 
in agreement with our earlier result [24]. Typical phase diagram of SL states in the (x, 5) 
plane with V = 2 is also shown in fig. (6). The phase diagram contains one to seven SL 
states. Furthermore, if we compare the results for 5 = with 5 ^ 0, we note that maximum 
number of possible SL states increases as nondegeneracy is introduced. 

In passing it is worth to mention that there are N/2 stable SL states in a region if the 
region in the phase plane contain N number of SL states and N is even. On the other hand 
(N + l)/2 states are stable if N is odd. The stability of a SL state is connected with the 
variation of the energy of the state as a function of x- If the energy of the state increases 
with the increase of x, the state is stable otherwise unstable. A simpler analysis of stability 
of a state and its connection with its energy will be presented in the next section. 

III. FULLY NON-LINEAR CHAIN WITH ALTERNATIVE NONLINEAR 

STRENGTHS 

Here we consider a fully non-linear 1-dimensional system where the alternative sites are 
of different strengths. The Hamiltonian of the system is given by 

1 OO 1 oo 

H =—7 1 £ Xn\CX +2 + - £ [C* n C n+1 + h.c] (13) 

' n=— oo n=— oo 

where , 

X2n = Xi, - oo < n < OO 

X2n+1 = X2, ~ oo < n < oo. (14) 

Hopping matrix elements between neighboring sites are taken to be unity. To obtain SL 
states we consider C n = (p n e~ lEt and n = 0O 7 ?'™' • The expression for i] is obtained by 
taking \n\ — > oo and is given by r] = ±[\E\ - s/E 1 - 4] [21]. Here we have considered 
monomeric ansatz because the system is symmetric about the 0-th site. After going through 
the same procedure as earlier and using the normalization condition , 0q = (1 — i] 2 )/2 , we 
obtain the reduced Hamiltonian of the dynamical system given by 
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H, 



1 



eff 



1 — 77 2 
1 + ?7 2 



<t+2 



1+ „a+2 l_„^+2 

x- — + 



+ 



2 V 



1+V 2 



(15) 



where x an d ^ are defined as x — (Xi + X2)/2 and 5 = (xi — X2)/2. Setting dH e ff/dr) = 
we get 

1 (rf - l)(?f+ 4 + l) 



X 



rp-l f 1+r; 2 \ 2 _ r(rT+l 



)(>? CT + 4 -l) 



(1 — ^+2)2 



(16) 



(1 +t; <t+2)2 

Eq.(16) can be used to analyze number of possible states for different values of x, 8 an d a. 
For a = and 5 = from eq.(16) we see that to get a SL state x needs to be infinite. It, 
therefore, means that no SL states can be obtained. This is true for a = and 5 = 
because the system reduces to a perfect linear system. This is also true for a = and 5^0. 
Analysis of SL states can be done for arbitrary a but we consider the case of a = 2 for the 
reason mentioned earlier. For a = 2, eq.(16) reduces to 
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4\2 



?7 (l + ?7 6 )(l + ?7 4 ) 



(17) 



X (1 + ^ 2 ) 3 [(1 - ?? 2 )(1 + V 4 ) 2 ~ <fy(l " ^ 2 )(1 - V 6 )} 
Eq.(17) directly tells that for 5 = there will be always one SL state and this is consistent 
with our earlier result [24]. Using eq.(17) the phase diagram of SL states in (x, S) plane 
is obtained and shown in fig. (7). As usual numbers in different regions indicate number of 
SL states in those regions. Here also maximum number of SL states increases compared 
to a perfect nonlinear chain. Furthermore the phase diagram in this case is quite rich (see 
fig-(7)). 

The stability of the SL states in this system can be understood from a simpler graphical 
analysis. For this purpose we look at the Hamiltonian given in eq.(15). The Hamiltonian 
is function of one dynamical variable, namely, i] where x an d 5 are parameters. For fixed 
values of x an d 5, rj and dH e ff/drj = f(rj) can be treated as generalized coordinate and 
generalized momentum respectively for the reduced dynamical system. Solutions of the 
equation f(r)) = will give the fixed points for fixed value of x an d 5 and number of fixed 
points with rj G [0, 1] will be the number of SL states. We now plot f(rj) as a function of rj 
for x — 0.35 and 5 = 2.5 as shown in fig. (8). If f{rj) > 0, then the flow of the dynamical 
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variable will be along the positive direction and on the other hand if f{rj) < the flow will 
be along the negative direction. These flows are indicated by arrows in the figure and the 
fixed points are denoted by A, B and C respectively. In the neighborhood of A the flow is 
always towards A and same for C. Therefore, A and C are stable fixed points. On the other 
hand in the neighborhood of B the flow is away from B. Hence, this is an unstable fixed 
point. Same thing happens for all (x, 5) in the three state region in fig.(7). We therefore 
notice that among three SL states (fixed points) two are stable and the other one is unstable. 
To get a connection of the stability of the states with the energy and x, we plot the energy 
of the states as a function of x in the neighborhood of x — 0.35 for fixed value of 5 = 2.5 
in fig. (9). Energies of the SL states arising from the fixed points A, B and C are denoted 
by Ax, Bi and C\ respectively in fig. (9). The figure clearly shows that in the neighborhood 
of x — 0.35, the energy of two stable SL states increases with x an d that of the unstable 
state decreases. We, therefore, conclude that the SL state is stable if the energy of the state 
increases with x otherwise unstable. 



IV. SUMMARY 

We have found the possible number of SL states due to the presence of a quadratic 
nondegenerate nonlinear dimeric impurity in one dimensional systems. Phase diagram of 
SL states in (x, 5) plane is presented for all the systems with a nondegenerate dimeric 
impurity. Maximum number of SL states is found to increase due to the introduction of 
nondegeneracy. The dimeric ansatz has been used in this case and this reproduces known 
results for special cases. A full phase diagram of SL states for a system of one dimensional 
chain with alternative sites having different nonlinear strength is presented. A monomeric 
ansatz has been introduced in this case and found to be consistent for special cases. Here 
also maximum number of SL states increases compared to that in a perfect nonlinear system. 
A stability analysis for the SL states in the fully nonlinear system and a connection of the 
stability with the variation of the energy of the state as a function of the nonlinear parameter 
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are presented. 



10 



REFERENCES 

[1] Holstein T D 1959 Ann. Phys. 8 325, 8 343; Turkevich L A and Holstein T D 1987 
Phys. Rev. B 35 7474 

[2] Eilbeck J C, Lomdahl P S and Scott A C 1985 Physica D16 318 

[3] Davydov's Soliton Revisited: Self- Trapping of Vibrational Energy in Protein Vol. 243 
of NATO ASI Series, Series B: Physics edited by Christiansen P L and Scott A C 
(Plenum, New York, 1991). 

[4] Emin D 1982 Phys. Today 35(6) 34 

[5] Toyozawa Y 1983 Organic Molecular Aggregates edited by Reineker P, Haken H and 
Wolf H C ( Springer- Verlag, Berlin) 

[6] Campbell D K, Bishop A R and Fesser K 1982 Phys. Rev. B 26 6862 

[7] Chen D, Molina M I and Tsironis G P 1993 J. Condens. Matter 5 8689 

[8] Wan Y and Soukoulis C M 1989 Phys. Rev. B 40 12264 

[9] Molina M I and Tsironis G P 1994 Phys. Rev. Lett. 73 464 

[10] Kenkre V M, Tsironis G P and Campbell D K 1987 Nonlinearity in Condensed Matter 
ed. Bishop A R, Campbell D K, Kumar P and Trullinger S E (Springer- Verlag) 

[11] Cohen M H, Economou E N and Soukoulis C M 1983 Phys. Rev. Lett 51 1202 

[12] Chen W and Mills D L 1987 Phys. Rev. Lett. 58 160 

[13] Abrikosov A A, Buzdin A I and Kulic M L 1989 Supercond. Sci. Technol. 1 260 
[14] Claude C, Kivshar Y S, Kluth O and Spatschek K H 1993 Phys. Rev. B 47 14228 
[15] Dauxois T, Peyrard M and Bishop A R 1993 Phys. Rev. E 47 684 
[16] Cai D, Bishop A R and Gronbech- Jensen N 1994 Phys. Rev. Lett. 72 591 

11 



[17] Tsironis G P, Molina M I and Hennig D 1994 Phys. Rev. E 50 2365 ; Molina M I and 
Tsironis G P 1993 Phys. Rev. B 47 15330 

[18] Yiu Y Y, Ng K M and Hui P M 1995 Phys. Lett. A 200 325 ; 1995 Solid State Commun. 
95 801 

[19] Hui P M, Woo Y F and Deng W 1996 J. Condens. Matter 8 2011 

[20] Gupta B C and Kundu K 1997 Phys. Rev. B 55 894 

[21] Malomed B and Weinstein M I 1996 Phys. Lett. A 220 91 

[22] Aceves A B, De Angelis C, Peschel T, Muschall R, Lederer F, Trillo S and Wabnitz S 
1996 Phys. Rev. E 53 1172 

[23] Gupta B C and Kundu K 1997 Phys. Rev. B 55 11033 

[24] Gupta B C and Kundu K (unpublished); Kundu K and Gupta B C (unpublished) 
[25] Economou E N 1983 Green's Function in Quantum Physics (Springer- Verlag, Berlin) 



12 



FIGURES 

FIG. 1. The phase diagram of SL states in S) plane for an one dimensional system with a 
nondegenerate linear dimeric impurity in the middle of the system. The (%, S) plane is divided into 
many regions by the critical lines. The region indicated by 2LA implied the possibility of two SL 
states in that region and both of them lie above the band of the linear host. Similarly 1LA, 1LB 
and 2LB implies one SL state above the band, one SL state below the band and two SL states 
below the band respectively. 

FIG. 2. The phase diagram of SL states in 5) plane for a Cayley tree with a nondegenerate 
linear dimeric impurity embedded in the middle of the system. The connectivity (K) of the Cayley 
tree is 4. The region indicated by OL implies that there is no possible SL state in that region. 
Others are same as in fig. (1). 

FIG. 3. The phase diagram of SL states in (x, 5) plane for one dimensional system with a 
bond defect (V = y/2) between sites and 1 as well as a nondegenerate linear dimeric impurity 
occupying sites and 1. Other things are same as fig. (1) and fig. (2). 

FIG. 4. The phase diagram of SL states in (x, 5) plane for an one dimensional system with a 
nondegenerate quadratic nonlinear dimeric impurity in the middle of the system. The (%, 5) plane 
is divided into many regions by the critical lines. The number in a region indicate the number of 
possible SL states in the region. 

FIG. 5. The phase diagram of SL states in (x, S) plane for a Cayley tree with a nondegenerate 
quadratic nonlinear dimeric impurity embedded in the middle of the system. Here K = 4. The 
number in a region indicates the number of possible SL states in the region. There are two small 
regions near the origin containing two SL states. 
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FIG. 6. The phase diagram of SL states in (%, 5) plane for one dimensional system with a 
bond defect between sites and 1 as well as a nondegenerate quadratic nonlinear dimeric impurity 
occupying sites and 1. The bond defect, V is taken to be 2. The number in a region indicates 
the number of possible SL states in the region. 

FIG. 7. Phase diagram of SL states in (x, S) plane for a fully nonlinear one dimensional system 
where alternative sites are of different nonlinear strengths. The number in a region indicates the 
number of possible SL states in the region. 

FIG. 8. f{rj) is plotted as a function of r/. Points indicated by A, B and C are fixed points. 
Here % and 5 are taken to be 0.35 and 2.5 respectively. 

FIG. 9. Energy of SL states are plotted as a function of %. Here 5 = 2.5. A\, B\ and C\ are 
the points corresponding to the points A, B and C ( in fig. (8) ) respectively. The dotted line is 
X = 0.35. 
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